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Slotted Test Section Dynaalcs of the NTF 


1. Introduction . 

In the previous reports*, eoneerninig the aodelling of the 
NT? dynamics, a model of the wind tunnel system proper toget- 
her with a simplified control system has been presented. As 
a continuation of this effort a dynamic model of the test 
section and plenum has been developed. In arriving at an 
acceptable representation the following objectives were set 4 

a. The model should serve as a means for the determining 
the dynamic performance of the NT? and the subs^qunht design 
of optimal control configurations. 

b. The model should be of the one dimensional flow type 
so as to facillitate computer simulation. 

c. The model should allow for the similation of the effects 
of test section blockage, boundary layer losses, slot flow 
required for supersonic flow generation, reentry flap pos- 
ition, and auxilliary plenum suction. 

d. The model should be capable of validation by compar- 
ison to previous tests made on the eight foot tunnel at 
Langley. 

e. The model should be capable of being updated as oper- 
ating data from the NT? becomes available. 

f. The model should lead heuristically to simpler models 
so that the principles of operation of this part of the 
tunnel are more clearly understood. 

1.2 The model presented in this report meets most of the 
objectives. The shortcomings of the current representation 
are; 


•References ^-8 
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a. There is no representation of reentry flap position. 

b. The boundary layer coefficients Cxa are not known to 
a sufficient degree of precision. 

c. The test section divergence is not represented. 

d. There is no simplified representation available so as 
to better understand the test section dynamics. 

e. The actual value of Ar is not known. Ar is an equiva- 
lent area and may not be equal to the geometric area. Ar 
may even vary with mach number. 

f. The assumption of a constant loss coefficient for 
the diffuser may not be a good one. 

2. Basic Configuration . 

The basic configuration of the slotted test section 
and plenum is illustrated in figure 1. The following equa- 
tions are used. 

2.1 The flow entering the test section is given by: 


144*Atest»Pstag 


Qthrt = 


< 


( 1 + .2-Mach^) ^ 


•Mach / 1.4* 

5 yTZ 


e 

Tstag 


if Mach^l 


144- Atest* Pstag 


/ 1.4- 

Y R»Ts 


8 

Tstag 


if Mach > 1 


where : 


Atest is the test section area, Atest= 66.77 ft* 
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R is the imlversal gas constant, R« 55*1 

2 

g is the gravity constant, g= 52.2 ft/sec 

Pstag is the stagnation pressure (Ib/in ) 

Tstag is the stagnation temperature C*R) 

Mach is the test section mach number 

2.2 The mach number is given by: 

Mach = 

where : 

2 

Pstest is the static test section pressure (Ib/in ) 

2.5 The flow (Ibs/sec) entering the plenum from the test 
section via the slots is given by: 

Qslot = Qthrt* Crm + (Qthrt-Qtest) + Qdm 
where : 

^b = Qthrt»Cxm is the component due to boundary layer. 
Cxm is a function of mach number (see section 3«) 

Qt = (Qthrt-Qtest) is the component due to supersonic 
flow generation. 

Qdm = Qthrt-Qdmx is the component due to model blockage. 
Qdmx is an input variable. 

2.4 The flow exiting the test section is given by: 

Qexit = Qthrt - Qslot 
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2,5 Qtest is the test section flow and is given by: 


Qtest = 


144 «Atest* Pstag .Mach 
(1 + ,2-Mach^) ^ 



1.4.g 


R'Tstag 


2.6 The flow leaving the plenum and flowing into the diff- 
user via the reentry flaps is given by: 


Qflap = Qr - Qexit 
where : 

Qr is the total flow into the diffuser. 

2.7 The total weight of gas in the plenum is given by: 


Wpl = 



(Qslot - Qflap) dt 


2.8 The static pressure in the plenxim is given by: 


Wpl» R»Tpl 

Ppl = 

144. Vpl 

where : 

Tpl is the plen\im static temperature. It is assumed 

that the total and static temperature in the plenum 
are equal to each other as well as to the test 
section total temperature. 

z 

Vpl is the plenum volume, Vpl = 56,000 ft 
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2.9 The test section static pressure is assumed equal to 
the plenum static pressure. 

Pstest = Ppl 


2.10 The static temperature in the test section is given 
by: 


Ttest 


Tstest 


(1 +.2 Mach^) 


2.11 The velocity in the test section is given by: 


Utest = Mach Vi ,4- g .R .Tstest 


2.12 The elector action of the diffuser is represented by: 

144»Ar. Pstest + Qexit ►Utest/g = 144*Ar.Pr (1 +1.4-»Mr^) 
where : 

The balance is taken across an equivalent area Ar. The 

exact value of this area is not known. A nominal value of 

2 

Ar= 1.16 Atest = 77.4-5 ft was chosen. The value of Ar was 
varied during the simulation so as to obtain a measure of 
the sensitivity of the system to the value of Ar. This is 
the familiar momentum equation applied across a line. The 
equation balances the sum of the potential and kinetic 
forces on either side of a line. On the left face of Ar 
the potential force is 144- .Ar -Pstest , The kinetic force 
is given by Qexit *'Jtest/g, since we may assume that the 
kinetic energy associated with the flow Qflap is negligible. 
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The energy on the right face of Ar is obtained by using 
the identity. 

144»Ar‘Pr.g 1.4 Mr^ = Qr»Ur 

2.15 The remainder of the tunnel » up to the fan, is repre- 
sented by a sequence of momentum, energy, and continuity 
equations. Ihe first equation encountered is: 

= [I'W-Ag.g (Pr - Pj^Q -IqgQgOg) + Qrur - /XL, 

where : 

The equation is written across a section of the diffuser 

of length XLg, XL^ =30ft, and an equivalent area of A^, 

2 

Ag =90ft . The state of the gas entering the section is 
(Pr,Qr,Ur). The state of the gas leaving is (^10*^10^10^* 

Qr and are assumed equal. This assumption is valid since 
it is equivalent to saying that the flow entering a volume 
of short length is equal to the average flow in that volume. 
Xq^ is the pV loss coefficient in this section. 

Xq^ = .1111 10“^ sec^/Cin^ft) 

2.14 The equations of 2.1 to 2.13 are solved simultaneously 
in an iterative manner. Care must be taken in the calculations 
so as to remain on the proper branch of the solution. This 
will assure that choking and the shock equation are properly 
represented. The procedure used is outlined in Appendix A. 

In general the procedure proceeds as follows: 
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Let 


Pm * 144 ’Pstest «Ar + Qexit *Utest/g 

and 

d = Qr^R Ttr /(1.4«Pm^. g) 

whene Ttr is the total temperature on the right face of Ar. 
If there is no loss of heat in the plenum Ttr= Ttest'i . 

Mr is given by the quadradic expression. 

(5-^9 d)*Mr^ + (25- 70d)*Mr^ - 25 d = 0. 

The remaining variables are given by: 

Pr = Pm/(144. Ar(l +1.4Mr^)) 

Tr = Ttr/(1 +.2Mr^) 

Ur = 49.84 ‘Mr^/fF 

If Mr^l, choking occurs and we set Mr equal to 1. 

Qr is then given by: 

Qr = Fm g/(24.R.Ttr) 

. • 

Qg is set equal to Qr, and Qg is constrained so that Qg$0. 

5. The distribution of losses around the tunnel . 

In order to properly model the tunnel, the losses around 
the tunnel must be properly represented. From the report by 
D.M. Rao (11), it may be inferred that when operating at 
14.7 psia, 500^R the tunnel drops will be disillbuted as 
follows; 

station 1 to 4 13.91% 

station 4 to 6 3.21% 
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station 6 to 8 11.37% 
station 8 to 10 60.86% 
station 10 to 12 5«^% 
station 12 to 14 5»19% 


From the report by B. Gloss (2) the fan ratio as a function 
of mach number is given as; 


Mach No . 

:2 

.5 

.4 

.5 

.6 

.7 

.8 

.9 

1.0 

1.1 

1.2 


Fan Ratio 

1.0120 

1.0180 

1.0243 

1.0350 

1.0450 

1.0550 

1.0700 

1.0880 

1.1120 

1.1420 

1.1850 


3.1 Using the above data the losses around the tunnel may 
be established if we further assxune: 

2 

a. The losses arepv losses, with constant loss 
coefficients. 

b. That the slot flow due to boundary layer at 
Mach 1 is 5% of the flow at the throat. 

c. That the slot flow is a function of mach number. 

5.2 The coefficients Xq^, i=3, 5, 7,9, 11,13 are defined by; 

DPi = Xq^.Q..U. 

where ; 

DP^ is the pressure loss in the ith section. 

Xq^ is the loss coefficient in the ith section. 
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is the flow in the ith section, 
is the velocity in the ith section. 

The coefficients Cxm^ are defined by: 

Qb= Cxm^- Qthrt 

where : 

The subscript i denotes the mach number at which the 
computation is made. 

Qb is the slot flow due to boundary layer. 

Cxm^ is the slot flow coefficient. 

3.5 The coefficients Xq and Cxm are determined as follows: 

a. Establish the tunnel simulation at 14.7 psi, 500 
Mach 1.0. Fix Cxm at .05 and vary all of the Xq^ in an 
iterative fashion until the fan ratio is 1.112 and the 

losses are distributed as indicated in 5.0. 

b. Establish the mach number at .1 to 1.2 in steps of 
.1 and for each mach number vary Cxm^ until the proper 
fan ratio is obtained. 

If this is done the following values are obtained. 

Xq^ = .8055 lO"*^ 

Xq^ = .3491 lO"*^ 

Xq^ = .1934 10 "^ 

Xq^ = .1111 10"^ 

Xq^^ = .1791 lO"*^ 

Xqj^^ = .2013 lO”*^ 
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Mach No* 

Cxm 

.1 

.4501 

*2 

.1612 

.5 

.0807 

*4 

.0445 

.5 

.0300 

*6 

*0232 

.7 

*0206 

*8 

.0216 

.9 

*0243 

1*0 

*0300 

1*1 

.0299 

1*2 

.0216 


3*4 The data given in 5*5 is suspect 
of Cxm at the low mach numbers* This 


due to the large values 
is more than likely due 


to the l^biccuracies of the assumed fan ratios. In order to 


correct for this error we introduce the energy ratio* 


Qthrt Utest^ 

CE = 

2(550)g HP 

where HP is the horsepower needed to drive the tunnel. 

The determination of the coefficients Cxm is redetermined 
Below mach numbers of .6 the fan ratios as given in 3*0 are 
ignored. Instead the values of Cxm are chosen so that the 
value of CE is held constant and equal to that at Mach =*6* 
This results in coefficients which appear to be more realis- 
tic, It is these values that are to be used to represent the 
test section* The results of this calculation is tabulated 
on the folowing page* 
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Mach No. 

Cxm 

CE 

Or 

Crl 

.1 

.0137 

4.70 

1.0012 

1.0080 

.2 

.0144 

4.70 

1.0048 

1.0120 

.5 

.0155 

4.70 

1.0107 

1.0180 

.4 

.0169 

4.70 

1.0189 

1.0245 

.5 

.0210 

4.70 

1.0299 

1.0350 

c6 

.0232 

4.70 

1.0430 

1.0430 

.7-. 

.0206 

4.73 

1.0550 

1.0550 

.8 

.0216 

4.96 

1.0700 

1.0700 

.9 

.0245 

4.89 

1.0880 

1.0880 

1.0 

.0300 

4.65 

1.1120 

1.1120 

1.1 

.0299 

4.59 

1.1420 

1.1420 

1.2 

.0215 

4.01 

1.1850 

1.1850 


Where Or is the resulting fan ratio and Crl is the fan ratio 
given in 3.0. 
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The X intercept of .35 corresponds to a dead time of .55 
seconds. This agrees with the acoustic propagation time of 


5^3 ft/(863 ft/sec) » .397 sec 

It may be shown that the slope of the line is the reciprocal 
of the time constant. This yields a time constant of 
1/.411* 2.43 sec. If the time constant is determined in the 
usual fashion, by defining the time constant as the time 
required to reach .632 of the final change, a value of 2.8 
seconds is obtained. Note that 2,43 + .35= 2, 76* 

Prom this we conclude that one can model the test section 
response by: 


AM = 


(TC‘S + 1) 


AG 


where : 

G is the change in guide vane angle. 

M is the corresponding change in Mach number. 

Tc is the time constant 
Td is the dead time, 
k is the gain. 

S is the Laplace transform operator. 


In designing the mach number control system, the above 
parameters prove useful. In the next section we examine how 
these parameters vary with the ope'-ating point, It is seen 
that they are a function of pressure, mach number, and Ar. 
In addition a nonlinear affect is observed in that the time 
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constant is dependent on the step size. Further for small 
value:- of AT choking occurs at supersonic operation. This 
choking occurs at che boundary of Ar amd causes velocity 
limitting which in turn slows down tbe response. 

h.2 The variation of gain and time constant. 

On figures 2a and 2b the gain and time constant of the 
test section is plotted as a function of stagnation pressure 
and aach number. This data was obtained by runs made on the 
STP simulation. With the pressure and temperature control 
active a 2 degree step change in guide vane angle was made. 
The time constant was determined as the time required for 
the mt-ch number to make .652 of the total change. Since 
we have already shown that the system is indeed a first 
order, this calculation is valid. However the time constants 
shown are the true time constant plus the system dead time. 
Since the distance between the test section and the fan is 
approximately 5^5 ft, the dead time varies between .507 and 
.487 seconds for the temperature range of 500 to 200 *E. 

To obtain this data runs were made at pressures of 30,60,90, 
and 120 psi; temperatures of 200,500,400, and 500 *R; and 
r.ach numbers of .2, .5, .4, . 5 , .6, .7, .8, .9,1.0, and 1.1. 

There was no change in either gain or time constant with 
temperature. Since it is known that the dead time does 
vary with tem;erature, it must be true that the actual time 
constant also varies with temperature in such a manner as to 
maintain the apparent time constant (qs plotted in fig 2a) 
fixed. The gain plotted on fig. 2b was obtained by taking 
the ratio of mach number change and guide angle change and 
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multiplying by 100 


4.3 The time constant variation due to Ar . 

It was observed that the value of Ar had a strong effect 
on the time constant. The curves of figure 4 show the vari- 
ation of time constant with Ar. It is noted that for low 
values of mach number that the value of Ar has little if 
any effect on the time constant. The time constant for tran- 
sonic operation is strongly dependent on Ar. These curves 
clearly show the need for obtaining be^er information as 
to the value of Ar. 

4.4 The time constant due to blockage . 

On figure 5, the time constant due to a 2% blockade is 
plotted. Comparing this with fiure 2a, it is seen that the 
the time constant for blockage and guide vane deflection is 
almost the same. That this is so is clearly seen by examin- 
ing figures 5 and 6. From these we conclude that the response 
due to blockage differs from that due to guide vane deflect- 
ion only in chat blockage does not have any dead time asso- 
ciated with it. Note that figure 5 has an x intercept of 
.35 seconds, while the line of figure 6 goes through the 
origin. 

4.3 Comparison with the eight foot tests . 

In order to assist in the validation of the NTP simula- 
tion transient tests were made on the eight foot tunnel at 
Langley. These tests were made in order to determine the 
response of the plenum-test section to test section block- 
age. The findings of these tests are documented in ref (10). 
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The table below whioh appears on page 11 of the aboTe report 
suimarizes the pertinent results. 


Test 

Run 

Nominal 
Kach No. 

SEL (#44) 

Average Time 
FM(IOIT) 

Constant, sec 
Eqs. (11) to (15) 

19 

1.05 

7.7 

9.2 

13.5 

20 

1.20 

5.7 

5.2 

3.5 

21 

.95 

7.1 

6.7 

4.5 

22 

.80 

5.5 

5.5 

4.6 


In order to cosipare the present NTP simiilation with 

these tests, the NTF simulation was altered so that it would 

exhibit the properties of the eight foot tunnel. This was 

accomplished by changing volumes and lengths in the tunnel 

proper. This change is described in ref (8). The plenum 

volume was increi^sed from 56,000 to 72,000 ft^.The value 

2 

of Ar was left at 77.^5ft , and a linear fan representa- 
tion was utilized. The inbleed and outbleed valves were 
locked and a 2% step change in blockage was made in the 
test section. . Pour runs were made : run 2219, Tstag=500, 
Pstag=l4,7i, Mach= 1.05; run 2220, Tstag=500, Pstag=14.71, 
Mach=1.2; run 2221, Tstag=500, Pstag=14.71, Mach= .95; 
and run 2222, Tstag=500, Pstag= lh.71, Mach= .8. These 
correspond to runs 19,20,21, and 22 respectively. A 
typical response is shown on figure 8, 

An analysis of the output of figure 8, is made by the 
use of figure 7» From this figure we conclude; 
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a. The response of nach nuaber is of the first order. For 
run 2219* with an Ar of 77.^5it^, the time constant is 
3.95 sec, there is no system dead time. 

b. The response of the test section pressure is of the first 
order. For run 2219, with an Ar of 77. the time constant 
is 3.93, there is no system dead time. 

c. The response of stagnation pressure, is of the third 

2 

order. For run 2219, with an Ar of 77.^5ft , there is a 
time constant of 3.69 seconds, with an additional underdamped 
second order of approximately 3 hertz. There is a system 
dead time of .33 seconds. 


The table below compares the eight foot tests with 
the simulation results. 


Run 

8 ft tests 

simulation 
Ar=70 Ar»72 

II 

Ar«77.45 

19 

7.7,9.2,13.5 

C 

6.5 

5.3 

4.0 

20 

5. 7, 5. 2, 3. 5 

C 

C 

19. 

2.9 

21 

7.1,6. 7,^. 5 

6.2 

5.0 

5.6 

3.9 

22 

5.5, 5.5, ^.6 

^.9, 

^.6 

4.4 

4.7 


The listed values are- measured time constants in sec- 
onds. A C indicates that the simulation failed to run due to 
choking at Ar. 

Examining the above, it may be concluded that the sim- 
ulation concurs with the test. In order to obtain better 
correlation or detect differences between the simulation 
and the tests, it will be necessary to obtain mere consis- 
tent test data and establish more precise values of Ar. 
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5.0 Mach Niaiber Contro l 

Proa the low band-pass of the aach number response) it 
is apparent that the guide vane actuator wl!.l be the limitt- 
ing factor on the response of the aach number control system. 

The process time constant is of the order of 5 seconds, 
while the actuator is a second order system of 6 radians 
with a damping factor of .8. Together these two yield a 

o 

phase angle lag of 130 at 3 radians, which limits the resp- 
onse to 3 radians per second. This is in agreement with the 
analysis presented in (6). 

The settling time, to within 19^, for a second order 
system with a natural frequency of 5 radians per second 
and a damping factor of .8, is obtained by solving. 


e“»Q(3)t 


01 


This yields a settling time of 2 seconds. 

This is not consistent with the system performance as 
exhibited by the NTF simulation. A typical simulation output 
is contained on figures 9a, 9b, and 9c. This run corresponds 
to a step change in mach number set point of .1, with the 
tunnel operating at 100 psi, 300 *!r and mach=.8. In this 
run a settling time of 20 seconds was required. This is 
far different from that expected for a properly tuned 
system. The discrepancy is due to croos coupling, the 
presence of reset in the control, improper tuning and 
perhaps some nonlinear affect unaccounted for. This dis- 
crepancy is to be investigated in a later study, meanwhile 
the target values of 2 to 5 seconds for settling times of 
the mach number control system is to be maintained 
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6,0 Conclusions and Re commend at ions . 

From the results presented in this report, it may be 
inferred that the proposed test section math model is realise- 
tic. It is therefore reasonable to proceed with the dynamic 
investigation of the NTF utilizing the model established in 
this report. 

The following action items have arisen as a result of 
this study: 

1. Obtain better data as to the distribution of losses 
in the tunnel. 

2. Establish more precise values of fan ratios as a 
function of mach munber. 

3. Obtain better data so as to establish the boundary 
layer losses. 

Obtain data as to the variable losses in the high 
speed diffuser. 

Determine more precise values for the parameter Ar. 
Establish any variation in Ar with mach number. 

6. Design an experiment on an operating tunnel to 
facillitate further validation of the simulation. 

7. Proceed with the design of a suboptimal mach nximber 
control system. 
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Appendix A. The solution of the diffuser-e.jector equation 


1, The momentum balance across the area Ar is given bj: 
144 Ar Pstest + Qexit Utest/g = 144 Ar Pr (1+1.4 Mr") 

T. The momentum equation between Pr and P^q is given by: 
Qg =^l -4*A^»g (Pr Ug) +Qr»Ur 

5. To solve these, in terms of Pstest, Mach, and we 

write : 

Ur 

Mr= 

•^1.4»R*g.Tr 

Qr*(5 +7 Mr^) 

Fm *5(5 +Mr^) 

where - 

Fm »144 .pstest «Ar + Qexit* Utest/g 


Qr 


Ar* ^ ^1.4»R »g .Tr 


Irtr^R 

^ l.^-B 


Letting 

d = 


Qr^* R .Ttr 
1.4 *Fm^. g 


the above equation may be written as: 


5*Mr^ (5 + Mr^) 
= d 

(5 + 7*Mr^)^ 


30 



or 


(5 -49 d) + (25 -70 d) - 25 d = 0 


4, The above quadratic may be shown to satisfy: 
case 1 . d ^ 5/^9 

2 2 

There is one real solution for Kr . Further Mr 4 1. 
case 2 , 5/^9 < d < 5/2^ 

2 

There are two real solutions for Mr . Moreover if x is a 

2 

solution for Mr , then y is the dual solution. With x and y 
satisfying: 

5 + y 5 + X 

X = ; and y = 

7y -1 7x - 1 

further if x^l, then y^l, and visa versa, 
case 3 » d 1^5/24 

2 

For d= 5/24, the solution is Mr =1. For d>5/24 there 
is no real solution. This is interperted as choking and the 
computation is so arranged to constrained $5/24, 


5 . The above system is solved by letting: 


b 


25 - 70 d 25 d 

. c = 

2 (5 -'+9 d) (5 -49 d) 


Then, if d^5/^9, case one holds and 



b + 



c 
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For 5A9 < d < 5/2^* case two holds and 
Mr^ = b - Vb^ - c 

For d^^/2^y the solution is degenerate. In this case v;e set 
Mr equal to 1, The value of Qr is obtained from: 

^r = Fm 



6,0 The ejector equation is seen to be consistent with 
the normal shock equation, if Qr= Qexit, For this case we 
get due to the symmetry between Mach and Mr, 

5 Mr^ (5 + Mr^) 5 Kach^ (5 + Mach^) 

= d = 

(5 +7 Kr^)^ (5 +7 Mach^)^ 


If Mach ^1, then d ^ 5/24 and either case one or case 
two holds. For case one Mr =Mach, since there is only one 
solution. For case two, again Mr= Mach, since the other 



is not physically realizable ( Mr can not be greater than one). 


If Mach^l, t;hen 5/49 ^d <5/24, and case two holds. 
There are now two possible solutions 

Mr = Mach, and Mr = +Mach^^/( 7 Mach^-1) 

which is the well known shock equation. 
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